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Abstract We consider a heavy, uniform, elastic beam rested on periodically distributed 
supports as a simplified model of a bridge. The supports are subjected to a partial destruction 
propagating as a failure wave along the beam. Three related models are examined and 
compared: (a) a uniform elastic beam on a distributed elastic foundation, (b) an elastic 
beam which mass is concentrated at a discrete set of points corresponding to the discrete set 
of the elastic supports and (c) a uniform elastic beam on a set of discrete elastic supports. 
Stiffness of the support is assumed to drop when the stress reaches a critical value. In the 
formulation, it is also assumed that, at the moment of the support damage, the value of the 
'added mass', which reflects the dynamic response of the support, is dropped too. Strong 
similarities in the behavior of the continuous and discrete-continuous models are detected. 
Three speed regimes, subsonic, intersonic and supersonic, where the failure wave is or is 
not accompanied by elastic waves excited by the moving jump in the support stiffness, 
are considered and related characteristic speeds are determined. With respect to these 
continuous and discrete-continuous models, the conditions are found for the failure wave 
to exists, to propagate uniformly or to accelerate. It is also found that such beam-related 
transition wave can propagate steadily only at the intersonic speeds. It is remarkable that 
the steady-state speed appears to decrease as the jump of the stiffness increases. 
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1 Introduction 



In a system like a uniform or a periodic waveguide, it may happen that a locahzed damage 
causes a failure wave. Mechanically such process is similar to the phase transition, and it can 
occur if the wave is accompanied by a permanent energy release sufficient to overcome the 
associated energy barrier. A well-known example is a falling row of dominoes. Many studies 
have been devoted to the related dynamic problems. Plane crushing waves were consid- 
ered in Galin and Cherepanov (1966), Grigoryan (1967), Slepyan (1968, 1977), Slepyan and 
Troyankina (1969), where, in particular, the uncertainty of a condition at the wavefront was 
discussed. It was then shown that uniqueness can be achieved in the framework of a struc- 
tured material model where the speed- dependent wave resistance to the transition can be 
determined (see Slepyan, 2002). To describe the related phenomena a higher-order-derivative 
formulation for an elastic continuum was used (see Truskinovsky (1994, 1997), Ngan and 
Truskinovsky (1999), and the references therein), and a discrete chain model (Slepyan and 
Troyankina (1984, 1988)). Waves in discrete bistable chains were then studied in Puglisi and 
Truskinovsky (2000), Slepyan (2000, 2001), Balk et al. (2001a,b), Cherkaev et al. (2005), 
Slepyan et al. (2005), Vainchtein and Kevrekidis (2012). Localized transition wave in a two- 
dimensional lattice model was considered by Slepyan and Ayzenberg-Stepanenko (2004), 
while non linear effects have been studied with large scale atomistic models by Buelher et 
al. (2003). 

With regard to a large-scale long-length construction, the failure wave may be supported 
by the gravity forces. In this connection we refer to papers by Baiant and Zhou (2002), 
Baiant and Verdure (2007), Baiant et al. (2008), which contain a comprehensive analysis 
of the collapse wave progress in the nine-eleven disaster. A bridge on pillars or a suspended 
bridge, an overpass, long conveyers are examples of the constructions where the failure wave 
may propagate taking energy from the gravity forces. 

In this paper, we examine some simplified models of such a construction considering the 
latter as a beam on a discrete supports and on a continuous elastic foundation, where the 
failure wave is that of a partial damage of the supports. Mechanically, these models differ 
from the above-mentioned ones by the existence of the subsonic and intersonic regimes of 
the failure wave propagation: no elastic wave is excited in the steady-state regime as far 
as it propagates in the subsonic speed range and there is such wave only behind the failure 
wavefront in the intersonic regime. Note that the elastic waves, if exist, create wave resistance 
to the failure wave. Possibilities of the wave radiation, in dependence on the failure wave 
speed, can be seen in the corresponding dispersion relations found for related structure in 
Brun et al. (2012). Dispersion diagrams to compare are plotted in Fig. [TDl 

The supported beam model also differs by the fact that the failure wave speed can vary 
in a wide range depending on the structure and damage parameters, and it can be very 
small compared with that in the above-mentioned bistable models. It is remarkable that 
the speed limit appears as low as the jump of the support stiffness is large. Three related 
models are examined: (a) a dynamic beam on a continuous elastic foundation, (b) a discrete 
set of masses rested on elastic supports and connected by massless beams and (c) a dynamic 
beam on the set of discrete elastic supports. Fig. [H Stiffness of the support is assumed to 
drop when the stress reaches a critical value. The failure wave is also considered under the 
condition that, at the moment of the support damage, the value of the 'added mass', which 
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reflects the dynamic response of tlie support, is dropped. We liave found tlie conditions for 
the failure wave to exists, to propagate uniformly or to accelerate. All speed regimes are 
considered, such where the failure wave is or is not accompanied by elastic waves excited by 
the jump in the stiffness and added mass. Limiting speeds are determined. We have shown 
that the effect of the change of the stiffness can essentially result in the failure wave speed 
limitations. 




Figure 1: Supported beam models. Undamaged and damaged parts are indicated with 
subscripts 1 and 2, respectively, (a) Beam on elastic foundation. The beam has mass density 
^1,2 = 1^12 + and bending stiffness D; the foundation has stiffness per unit length 

>^i,2 = x^i2/^- (b) Discrete set of masses = m^2^ + ^1,2 connected horizontally by 
massless beams of bending stiffness D and rested on elastic supports of stiffness x° 2 placed 
at distance a. (c) Dynamic beam of density mi2, added mass Mi^2 and bending stiffness D, 
on the set of discrete elastic supports of stiffness xj'g placed at distance a. The transverse 
displacement is denoted as w. 

In Fig. [2] we show an illustration of a bridge failed under a wave generated by an earth- 
quake. A periodic pattern is clearly visible in the collapsed structures and it corresponds to 
damaged foundations of the bridges. 
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Figure 2: View of the destroyed "Puente Viejo" over the Biobio river, that hnks Concepcion 
and San Pedro de la Paz, 500 km south of Santiago, Chile, after the 8.8 magnitude earthquake 
on February 27, 2010. (c) Nicolas Piwonka, National Geographic. 



2 Uniform continuous model: a beam on an elastic 
foundation 

We start with the simplest formulation of the structure as an infinite uniform elastic beam 
placed horizontally on a uniform elastic foundation as given in Fig. [1^. The foundation is 
considered as massless; however, the corresponding 'added mass' is assumed to be included 
in the beam per-unit-length mass, m. The foundation is initially stressed by the beam under 
gravity forces mg {g is the acceleration of gravity). We consider the failure wave as the 
propagation of partial damage of the foundation. Namely, we assume that the foundation 
stiffness drops at the wave front, t] = x — vt = 0, from its initial value, xi (ahead of the front) 
to X2 < Xi (behind the front), while the per-unit-length load on the foundation reaches a 
critical value, Qc at the front, rj = 0. Also, this problem is considered assuming that a part 
of the added mass disappears simultaneously with the jump discontinuity of the foundation 
stiffness. In this case, m = mi ahead of the front and m = m2 behind the front. Under the 
moving discontinuity, elastic waves can be excited propagating behind or/and ahead of the 
failure wave front; however, the static values of the beam vertical displacement are equal to 
rriig/ Xi (far ahead of the front) and to m2g/ >C2 (far behind the front). This allows us to 
make some conclusions based on energy considerations. 

2.1 Energy considerations 

Consider the loading diagram shown in Fig. [31 where gi 2, and Qc are the initial, the final and 
the critical load at damage acting on the support, respectively. The other values pointed in 
the figure are 



where wi^2 is the beam static displacement at 77 = ±00 and Wc is the critical displacement at 
?7 = -|-0. It can be seen that the work of the gravity forces, the static strain energy density 



X2 

q* = — Qc 

Xi 
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at ?7 = — oo and the energy loss due to the failure are 

1 11 

^ = + qi{wc - wi) + q2{w2 - Wc) , £2 = -^q2W2 , £* = ^{qc - q*)wc ■ (2) 

Since the energy of the beam itself in the steady-state regime is invariable, the energy excess 
per unit length is 



2 \X2 Xi y Xi 2xi \ Xi / 



(3) 



The energies are also shown in the load-displacement plane in Fig. [31 From this static value 
the energy radiated by elastic waves, if exist, must be deducted. In this general case, for the 
uniform propagating failure wave the energy balance is 

^0 - UM/v - 1) - U2{1 - C2/V) = , (4) 

where Ui^2 and C1.2 are the energy density and the group velocities of the waves at ±77 > 0. 
Note that f/i = in the sub- and intersonic regimes and f/2 = in the subsonic regime. 

In this problem, there is a range of the failure wave speed, Q < v < V21 where no elastic 
wave can be radiated in steady-state motion. The critical velocity ^2, indicated in Fig. HI 
is the resonant wave velocity corresponding to the lower dispersion curve. At this regime 
the phase and group velocities are equal. The value of V2 is given the next section. For 
this subsonic regime, we can conclude that the failure wave cannot propagate if £q < 0. In 
the opposite case, where there is a positive excess of the energy, the steady-state regime 
is impossible, and the failure wavefront move unsteadily spending the energy excess on the 
radiation of elastic waves, which arise under nonuniform motion. In the neutral case, Eq = 0, 
any value of the speed in the subsonic range satisfies the energy balance. However, in the 
latter case, the system parameters are connected by a relation following from the equation 
£^0 = 0. In particular, if no loss of the mass is assumed, this relation is 



Wc = W= Wi 



a/xi/^2- (5) 



In the subsonic regime, the wave propagates nonuniformly if the real critical displacement is 
less than this value, otherwise it cannot propagate. In the case £^0 > steady-state motion 
is possible at v > V2, solution of (jlj) as shown in the example at the end of Sect. [21 

Note that the last term in Eq. ([3]), the energy lost under the support damage, plays the 
role of the surface energy (or the effective surface energy) in fracture. The peculiarity of the 
considered problem is that the energy release under the support damage remains positive not 
only in the cases where there is no wave radiation, but also in a part of the intersonic speed 
range, and this occurs in the continuous model as well. In this latter case, the radiation 
decreases but not eliminates the total energy released under the 'phase transition'. 

2.2 Flexural waves, normalization and dispersion relations 

We based on the Bernoulli-Euler model for the beam and the Winkler model for the contin- 
uous foundation as in Fig. [T^. Thus, the dynamic equation is 

D +m — + Kw{x,t) =mg. (6) 
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Figure 3: Load q versus displacement w diagram for the beam on elastic foundation; qi and 
q2 are the gravitational loads for the undamaged and damaged structure of stiffness xi and 
X2, respectively, and Wi = qi/xi, W2 = q2/x2 are the corresponding static displacements. 
The critical load and displacement at damage are qc and Wc, respectively, whereas is the 
load corresponding to the displacement Wc of the damaged structure. A is the work of the 
gravitational loads, S2 is the static strain energy in the damaged structure at 77 — — 00 and 
is the energy loss at failure. The energy excess Sq is the difference between the area of 
the two triangles P^P^Pq and P1P2P3 in the {w, q) space. 

where w is the transversal displacement, respectively, k is the foundation stiffness and g is the 
acceleration of gravity. The foundation stiffness and possibly the density take different values 
ahead of and behind the wavefront. This continuous structure is assumed to be a limiting 
representation of the corresponding discrete-continuous system consisting of a continuous 
beam on a discrete periodic set of supports. Accordingly, the distributed mass density and 
support stiffness are assumed to be those averaged over the period, a 

m = m^ + M/a , x = x°/a , (7) 

where mP is the beam mass density, M is the added mass due to the supports and >iP is the 
discrete support stiffness (see also Sect. [3]). 

We consider the steady-state regime, assuming that the displacement depends only on 
T] = X — vt, where v is the failure wave speed. The equation of motion ([6]) becomes 

Dw{riY^ + mv^w"{r]) + >iw{rf) = mg (8) 
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with 



X = xi (?7 > 0) , X = X2 < xi (?7 < 0) , 

m = mi (?7 > 0) , m = m2 [r] < 0) . (9) 

In this equation, v is considered as the input parameter. Ones the steady-state solution is 
obtained, it must be used to satisfy the failure condition 

wiO) = w,. (10) 

If the latter equation can be satisfied it serves for the determination of the speed v; otherwise, 
the steady-state formulation failed as it was discussed in Sect. 12. ll 




Figure 4: Dispersion diagrams for the waves in uniform beams on elastic foundation of stiff- 
ness X. The dispersion curves are shown for the undamaged foundation {uJi{k) for stiffness 
X = xi), the damaged one {oJ2{k) for x = X2, with x = X2/X1 = 0.36) and the limiting 
case with x = (grey curve); no loss of mass is considered. Straight lines correspond to 
the different regimes of the failure wave speeds, v = uj/k= const. The group velocities 
C1.2 = doji 2/ dk coincide with the speed v at the resonant regimes, where v = V2 = \/2{ky^'^ 
and V = Vi = The (cu, /c)-points corresponding to the radiated waves are marked by small 
bullets; such resonant points are marked by larger bullets. The non-dimensional variables 
are used corresponding to the length and time units introduced in ffTTj) . 

We introduce the natural length and time units as 

i={DlH^fl\ T=^W/^l. (11) 
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After the normalization 

{x,r],w) = ^{x,f],w) , (t,w"^) = r(t,(I;"^) , {v, 01^2) = ^/r{v, 01^2) , 
k=~k/^, g={^/T^)~g, (£:o,^2,^*,Af/i,f/2) = Xie'(^o,^2,^*,At/i,f/2), (12) 

Eqs. (IE]) become 

w^^iv) + "^^ w"{rj) + w{rj) = g (77 > 0) , 
w^^ {ri) + mv'^ w'\ri) + kw{r]) = mg (rj < 0) , (13) 

where m = 7712/1711 and x = K2I y<\. Separating the initial static displacement of the undam- 
aged structure we implement the substitution w = w ^ g. Then, the above equation takes 
the form 

w^^(?7) +5^w"(?7) + w(?7) = (?7>0), 
w^^{ri) + mv'^ iv''{r]) + Hw{r]) = Q (?7 < 0) , (14) 

where Q = g [tti — k). 

In the following, the notations "tilde" and "bar" are omitted for ease of notation and 
we consider the case mi = m2. The displacements referred to in the following text are 
normalised dynamic perturbations with respect to the corresponding static values. 

The corresponding dispersion relations are 

Ui = ±VlT¥, (r/>0), U2 = ±VkT¥ (r/<0), (15) 

where k is the wave number. 

We consider three ranges of the velocity, v, corresponding to subsonic, intersonic and 
supersonic speeds. The boundaries separating these ranges correspond to equalities v = 
uj/k = ci^2 = (^^1,2/ (^k (ci^2 are the group velocities). We find that the three ranges are 
identified as follows: 

• subsonic range: Q < v < V2 = \P1{k)^I'^ , 

• intersonic range: V2 <'v <V\ = -\/2, 

• supersonic range: v > Vi. 

Dispersion diagrams ( 1T5|1 and sonic regimes are shown in Fig. HI Elastic wave radiation at 
the failure wave front 77 = corresponds to the intersection points of the dispersion relations 
( !T5|) with the propagating wave rays, u = v k, and the radiated wave propagates at 77 < if 
V > C2 and at ?7 > if < Ci. 

2.2.1 Subsonic regime 

Although, as follows from the energy considerations, the only 'neutral case' correspond to 
the steady-state regime in this speed range, we consider the subsonic regime in detail. We 
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show how the deformed shape of the beam in a vicinity of the transition point depends on 
the speed. Eqs. f lT3|) yield 

w{ri) = e"°^''(Ai cos Pit] + Bi sin (3ir]) (?7 > 0) , 
w{r]) = e"^''(A2 cos (32V + B2 sin /^sr/) + Q/k (r/ < 0) , 



as = ^V2v^-t;2, /32 = iy2Vx + y2^ (16) 

where the constants, Ai,...,B2, defined by the continuity conditions with respect to the 
displacement and its first three derivatives at 77 = 0, are 

Q Q 

^1 = ^=-^-=-, A2 



V^(l + y/k) ' k(l + y/k) 



V2 -v^{l + Vk) - 2V2Vk - 

-Dl = ^ I =~, ^ 2 : 



y/kV2Tv^(l - k) 



^ V27^^(l + v^)-2v^v/2^^ 
XV 2v X + v'^(l — k) 



The displacement at the transition point 

^(0) = r^, = - 1) g (18) 

Vx(l + Vx) VV ^2 / 

is independent of the speed as it should be in the subsonic speed range. It corresponds to the 
energy balance relation (|5]) (note that in (fT8|) we have separated the initial non-dimensional 
displacement equal to g). This regime can exist if the critical value of the displacement 
coincide with w(0) given by this relation. The plots of 'w{rf) and of w(0) as a function of 
X are presented in Fig. [5l In connection with energy considerations given in Sect. 12. 1[ we 
can see from Fig. [5b that, independently on the velocity v < a/2(x)^/^, the wave propagate 
nonuniformly if Xs/xi < {Wc/g + 1)~^ (in normalized coordinates) and it cannot propagate 
otherwise. In this regime the critical displacement for dynamic propagation coincides with 
the static displacement. 

2.2.2 Intersonic regime 

This regime is characterized by a sinusoidal wave propagating behind the transition point. 
Among two waves satisfying the equation we choose the wave with group velocity below the 
failure wave velocity (see Fig. H]), that is we keep the wave excited by the transition (the 
other propagating from minus infinity could correspond to a remote source). Thus 

w{ri) = e-"i''(Ai cos Pirj + sin /3ir/) (r/ > 0) , 

wijni) = A2 cos /32r] + B2 sin (327] + Q/ k (77 < 0) , 
1 



I32 = \jvy2-y/vy4-k, (19) 
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(a) (b) 

Figure 5: The beam on continuous elastic foundation. Subsonic speed regime (v < a/2x^/^): 
(a) Displacement profile w(?7), in the static case {v = 0) and at velocities v = 0.9, 1.15, for 
X = 0.5 (f2 = 1.189). (b) Displacement w/g at 77 = as a function of the stiffness ratio 
X = X2/X1. Dashed lines indicate the case shown in Fig. [5H. 



where the constants, Ai, ...,-82, are defined by the same continuity conditions as above for 
the subsonic case; they are 

v^-2k-Vv^-4k r- 2 ^'^ + 2>f - V^^^ - 4x 

A\ = ^ Q, B\ = —y2 — v^ , Q, 

2(1 -x)x -^'1 2v/2T^(l-x)x 



A2 = z — Q , B2 = -V2 -v^ ^ Q. (20) 

2 1-xx ^' ' v^(l-x)x ^ ' 



In particular, 



-2k- _ 

w{^) = ^-7 9- 21 

2x 



Since w(0) depends on f and this allows us to satisfy the condition (fTOl) . The displacement 
profiles wij]) for different values of the speed v are presented in Fig. [6ti. 

In Fig. [6b we also plot the value of w(0) as a function of velocity for x = 0.5. The 
steady-state propagation speed v corresponds to equality w(^) = Wc > ^ when X2/X1 < 
[wclg + 1)^^. We note that such a speed v depends on the safety factor wdg + 1 (expressed 
in term of non-dimensional values). 

As can be seen in Fig. [6)3 the displacement, ty(0), monotonically decreases as the speed, 
f , increases from the critical value, v = V2 = \/2(x)^/^, and becomes negative when the speed 
exceeds a value. These negative values correspond to a hypothetic case where the energy 
lost at the moment of the damage is less than its real minimal value, that is the critical load 
is less than the initial static load, qc < mg. 

The critical displacement w{0) versus v is shown in Fig. [7^ for different values of x, 
while the total speed bound w{0) = is given in the (x, f) space in Fig. [7)d. The steady 



10 



w/g w{0)/g 




(a) (b) 
Figure 6: The beam on an elastic foundation. Intersonic regime 

(a) Displacement profile ^(77) at different velocities, v = 1.19,1.25,1.40, for x = 0.5. (b) 
Displacement at 77 = as a function of velocity for the ratio x = 0.5 (1.189 < v < 1.414); 
the critical displacement Wc is attained at v = v*. 



propagation of the fault in the intersonic regime is not forbidden when w{0) > 0, which 
implies that the velocity v satisfies the following constraint 

< V < Vl + X. (22) 

It is also noted that for the intersonic velocities within the interval (-\/l + x, \/2) the failure 
wave propagation is not possible. 



2.2.3 Supersonic regime 

Finally, we address the case when the velocity of the transition front v > In this case, 
elastic waves are generated in both regions 77 > and 77 < 0. The displacement takes the 
form 

w{ri) = Ai cos /3i7/ + Bi sin /3i7/ (77 > 0) , 
w{r]) = A2 cos /32?7 + B2 sin l32r] + Q/x (77 < 0) , 

(3^ = ^vy2+^/vy4~l, (32 = \/ vy2 - ^v^i - k , (23) 

where the group velocity of the elastic wave ahead of the transition front is greater than the 
speed of the latter (see Fig. H]), and the constants defined by the continuity conditions are 

y/v"^ - 4 + y/v'^ - 4x X 

, + a s.=0. (24) 
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(a) (b) 

Figure 7: The beam on an elastic foundation. Intersonic regime (a/Sx^/^ < < y/2). (a) 
Displacement at 77 = as a function of velocity for the ratios x = 0.25, 0.5, 0.75, 0.9. (b) 
Total bound w{0) = at different values of the velocity v and the support stiffness contrast 
k. The grey area shows the intersonic regime, bounded by the dashed lines V2 = a/2(x)^/^ 
and Vi = \/2. Steady failure propagation is possible when w{0) > 0. 



Thus, 



v'^ — Jv'^ — 4x i — X , , 

w = , — , —9 < 0. 25 

Since the transition can occur when w{0) > this evidences that the wave cannot propagate 
at supersonic velocities. Thus, the absolute maximum of the failure wave speed is y/l + k. 

We conclude this section on the beam on elastic foundation with the results shown in 
Fig. [HI The displacement w{0) is plotted as a function of the failure velocity v for the 
stiffness ratios x = 0.35, 0.45, 0.75. It is possible to distinguish the subsonic regime (where 
w(0) is constant), the intersonic regime (where w{0) decreases with the increase of v) and 
the supersonic regime (where w{0) is always negative), delimited by the critical velocities 
V2 = a/2x^/^ and vi = \/2. At the critical displacement Wc/ g = 0.3 corresponds a maximum 
ratio X = 0.592 for the fault propagation. Therefore, propagation is not possible for x = 0.75, 
which is also in agreement with the fact that in (|3]) the energy excess So = —0.357(1 — 
k)g'^ /2 < 0. For x = 0.35, 0.45, where Sq is positive, steady-state propagation is given by 
the intersection points w{0) = Wc at v = 1.106, 1.164 (see Eq. (j2T])). respectively; two values 
which satisfy the constraint ( l22l) . It is trivial to verify that steady-state velocities satisfy 
energy balance (jl]) where 

C, = 0, U,^l-^(v'-v'V^P^), (26) 

4x^ V J V 



3 Discrete-continuous model 

A uniform elastic beam of the mass density mP is assumed to be placed on concentrated 
elastic supports at x = na,a > 0,?7, = 0,±1,... (see Fig. [It). We denote the stiffness, x°. 
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subsonic intersonic supersonic 




-0.6 L 



Figure 8: The beam on an elastic foundation. Displacement at the transition point is given 
in black continuous lines as a function of velocities in the subsonic, intersonic and supersonic 
regimes for Xi/x2 = 0.35, 0.45, 0.75. Grey line indicates the critical displacement wd g = 0.3 
and black circles correspond to the steady-state propagation of the fault at v* = 1.106, 1.164. 



of the intact and damaged supports by >c^ and Xg, respectively. Also the "added mass" M 
proportional to a is assumed to be placed at the junction points corresponding to positions 
of elastic supports. Thus, 

x° = x° , M = Ml , m° = nil (before the support damage, 77 > 0) , 
x° = X2 < x° , M = M2 , = ml (after the support damage, 77 < 0) . (27) 

The supported cross-sections have the coordinates x = na. The steady-state regime is 
considered in the sense that the displacement w at x = na is a function oi t] = na — vt, 
whereas the displacement of the beam between the supports is a function of two variables: 
w = w{x,r]). Thus, w{r]) = w{na,r]). 

The dynamic equations for the structure before the damage are 

D + ^1 = "^1^ {x^an)- 

Mi^^ + A^iji) - Q\ri) + Q~{ji) = Mig , 

M^{r]) - M'{r]) = {x = an), (28) 

where D is the bending stiffness, is the beam mass per unit length and g is the gravity 
acceleration; Q^{ri) and Ai^{ri) are the transverse force and the bending moment at the right 
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Figure 9: Inertial beam on elastic supports distributed at distance a. Panel I: configuration 
of the concentrated mass M. Panel II: auxiliary solution configuration for the beam. 



(superscript +) and at the left (superscript — ) of the junction point x = na, respectively 
(see Fig. [HI panel I). 



We introduce the same length and time units ^ = {D / Ki 



;V4 



\/m\l K\ as in ([TT 



with H\ = xj'/a, and the mass density unit mi = m\ + Mi/ a. In these terms, the normalized 
quantities denoted by tilde are 



(x, a, ?7, w) = ^{x, a, t], w), Mi^2 = miaMi^2 ■ 

, . d^w{x,ri) D~ ~ d^w{x,fi) 
Q{x, v) = -D = ^Q(x, v) , Q(x, v) = . 



M{x,r]) = D 



d'^w{x, 7]) 

dx'^ 



jM{x,r]), M{x,r]) = — 



(29) 



The equations of motion (1281) become 



d'^w{rj) 



ox* nil dt"^ mi 



Mi~—-^ + w{r^)-^ Q^{fi)-Q-{V) =Mig, (v) - (v) = {x = an) . (31) 

Recall that in the damaged region, ?7 < 0, the values x°,Mi, are replaced by X2 < x°,, 
M2 and For this case, equation (!30l) and the first equation in (!3T!) take the form 



d^w(x,ri) m^d'^wix.ri) mS , , . 

" H ^ ' " = —g {x ^ an) 

ox* mi ot^ mi 



M.^ + ^*M-ito-(f,)-Q-(f,) 



Xl 



M2g. 



(32) 
(33) 
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In the following the symbol "tilde" is omitted for convenience. 
We now separate the initial static state by the substitution 

w{x, rf) = g + ^ {x — anY{x — a{n + + w{x, rf) , (34) 

where the first term on the right-hand side represents the static displacement at x = na 
corresponding to the positions of the supporting pillars, whereas the second one stands for 
the relative displacement of the beam with respect to the supporting pillars and w is the 
normalised dynamic perturbation. Correspondingly, 

QHv) = ±!^^ + g±(,), mH^) = + A^^(^). (35) 

mi 2 nil 1^ 



We obtain that for r/ > 



d^w{x,ri) ^m1d'^w{x,ri) ^ 
dx'^ nil dt'^ ' 



M^{ri)- M-{r]) = Q, (36) 



x = an: Mi—^ + w{7])--[Q^{r])-Q-{rj)\=Q 



and for 77 < 



d^w(x,ri) d'^wixji]) — m? 

x^an: — — + — = g ; 

ox^ nil o't^ f^i 

x = an: M^^^^ ^ ^w{r^) -- {r^) - Q~ {t^)] = C , 
ot'^ Ki a ^ -' 

M-^{r])- M-{r]) = 0, 
C = g{\- X2/X1 + M2 - Ml) . (37) 

In the following text, the symbol "bar" are omitted for convenience of notation. We also 
restrict our attention to the case where the damage is concentrated at the pillars, positioned 
qX X = na, and = m2 = m°. 

3.1 Interaction of the neighboring supported cross-sections 

We now consider the interaction of the neighboring cross-sections placed at the elastic sup- 
port points, say, at a; = and at a; = a. For the steady-state regime the dynamic equation 
for the displacement of the beam placed between these cross-sections follows from fl36l) as 

d^wix.ri) ^wP d'^wix.ri) , , 

^ 4 " + ^4^ = 0- 38 

ox^ mi or]^ 

In terms of the Fourier transform on rj with parameter k we have 

aV(a:,fc) _ , witht;V = - lim(£ + ii;fc)2. (39) 

dx^ mi ^ ' ^ £^+0^ ^ ^ ^ 
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In the following, we will make use of an auxiliary solution of this equation corresponding to 
the boundary conditions (see Fig. [9l panel II) 



\0,k) = {w^{0,k)y = 0, w^{a,k) = wo, [w^ {a,k))' = (p^ . (40) 
It is 



W{x) 



= W{x)wo + $(x)0o 

(cosh Xa — cos Xa) (cosh Ax — cos Ax) — (sinh Xa + sin Xa) (sinh Ax — sin Ax) 

2(1 — cosh Xa cos Xa) 



^ (cosh Aa — cos Aa)(sinh Ax — sin Ax) — (sinh Aa — sin Aa)(cosh Ax — cos Ax) 

~ 2A(1 - cosh Aa cos Aa) ' ^ '* 



where A = (m^ /niiY^^y/kv — iO. The transverse shear force, Q = —w'", and the bending 
moment, Ai = w", follow from this as 

M^{0,ri) = M^oWo + M^o<Po, M^{a,ri) = M^aWo + M^a<Po, (42) 



with 



^ A^ (sinh Aa + sin Aa) ^ A^ (cosh Aa — cos Aa) 

""^ 1 — cosh Aa cos Aa ' 1 — cosh Aa cos Aa 

A^ (cosh Aa sin Aa + sinh Aa cos Aa) A^ sinh Aa sin Aa 

^wa Z r~r 7 ) ij>a 



1 — cosh Aa cos Aa 1 — cosh Aa cos Aa 

A^ (cosh Aa — cos Aa) A(sinh Aa — sin Aa) 



1 — cosh Aa cos Aa ' 1 — cosh Aa cos Aa ' 

A^ sinh Aa sin Aa A(cosh Aa sin Aa — sinh Aa cos Aa) 

1 — cosh Aa cos Aa 1 — cosh Aa cos Aa 

Note that for low speeds, w — )■ 0, the compliance coefficients tend to their static values, 
namely 

12 6 

QwO ~^ Qwa ~^ ! Q(hO ~^ Qsa ~^ 7^ ■, 

a-^ a^ 
6 4 

M^Q -Mu,a , M^a -2M^o - . (44) 
a^ a 

3.2 The Wiener-Hopf equation 

We now rearrange Eqs. (l36l) . (1371) making the Fourier transform and using expressions 
fH2E3D . We denote 

/•oo 

Jo 

{w.{k),c^.{k)}= I {^/;(r/),0(r/)}e''=''dr/, (45) 

J — oo 
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and deduce 



(1 - Miv^k'^)w+{k) + (X2/X1 - M2V^k^)w^{k) + -{Q^a - Qwo coska)w^{k) 

a 



2i C 

H — Q ^0 sin ka(j)^{k) = 

a + ik 



i-Mwo sin ka w^{k) + {Ad^a — A^</)0 cos ka)(j)^{k) = 



Excluding (t>^{k) we obtain 
Li{k)w+{k) + L2{k)w4k) 



C 



C{L,{k) - L^ik)) 



+ ik (0 + iA;)[l - xa/xi + (Mi - M2)(0 + ikvy] 



Li{k) = l + Mi{0 + ikvy 



L2{k) = ^ + M2(0 + ikvf + - 
Xi a 



{Qwa - Qwo COS A;a) - 
{Qwa - Qwo COS ka) 



Finally, the Wiener-Hopf type equation follows in the form 

C 

Lo{k)w+{k) +w_{k) 



Q^qMu^o sin^fcg 
M^a - A^</,oCos ka 

Q<f,oMyjo sin^ ka 
M^a - cos ka 



(0 + iA;)[l - X2/X1 + (Ml - M2)(0 + ikv)^] 

Lo{k) = L,{k)/L2{k). 



(46) 



(47) 



(48) 



3.3 Factorization 

To facilitate this action we assume that a small dissipation exists in proportion to the strain 
rate, that is, we correct the expression for the bending moment to be 



M. = T-^ h a- 



dtdx"^ 



(49) 



where a is a small time parameter of viscosity. The non-dimensional beam equation after 
the Fourier transform on rj becomes 



1^1 + ikva) 



d^w {x, k) 
dx^ 



v^k^w^{x,k) = 



(50) 



that results in the change of the parameter A 

A / \ y/kv — iO 



A 



— 



[l + ikvay/"^ \mi) (l + ifcwa)V4 
Now the product 

M 



0{\k\^'^) (t; > 0, -> 00) . (51) 



L{k) = j^Lo{k) ^ I (A; ^±00) 
has no zeroes on the real /c-axis (if X2 > 0). It means that its index is 0, i.e. 

IndL(/c) = ;:^[ArgL(oo) - ArgL(-oo)] = . 

ZTT 



(52) 



(53) 
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Besides, the integral of In L{k) over the real k-axis converges. This allows us to use the 
Cauchy type integral for the factorization of this function, that is to represent 



L(k) = lim LJk)Lik) , L±ik) = exp 
Eq. fHHl) can be represented in the form 



±- 



"27ri 



[±'^k > 0) . (54) 



L^{k)w4k) + = G{k) 



MiL_{k) 



M^^^^'^ - MAO 



G{k) 



MoC 



Mi(0 + iA;)[l - X2/X1 + (Ml - M2)(0 + ikv)^] ' 
First consider the case Mi = M2 = M. The latter function can be split by two terms 



(55) 



G{k) 



L+ik) 



L-{k). 



C1 + C2, Ci = ^[L+(fc)-L+(0)] 



C, 



9 



+ ik 



L+(0) - 



L4k) 



(56) 



where Ci,2 are regular in the upper and lower half-plane of k, respectively. The solution 
follows as 



w+{k) 



g L+(A;)-L+(0) 



ik 



LAk) 



W-{k) 



9 



Q + ik 

In particular, using the limiting relations 



-[L+(0)L_(fc)-l]. 



we find that 



with 



L±(0) 



w;(±0) = lim (^iA;)u;±(A;) , L±(±ioo) = 1 

fc— ^±ioo 



y;(+0) = w{-Q) = w{0) = g[L+{0) - 1] 
ArgL(/c) 




exp 



±- 



k 



dk 



m 



X2 



(57) 



(58) 



(59) 



(60) 



In the case M2 < Mi, the function G{k) in ( 155|1 can be split by the corresponding terms 
as follows 



Gi{k) = Go 



G{k) = Gi{k) + G2{k) 
L+(fc)-L+(0) _ 1 /L+(fc)-L+(/3) , L+{k) - L+{-/3) 
2 



C2{k)=G, 



+ ik 

L+(0) 1 



+ i{k-l3) 
LAP) 



+ 



+ i(A; + /3) 



+ iA; 2\0 + i{k-l3) + i(A; + /3) 
M2G 



G{k) 
L4k) 



Mi(l - X2/X1) 



1 / 1 - X2/ Xi 
^ ~ w V M1-M2 



(61) 



and the solution is obtained as 



l-L+(0)/L+(fc) 1 / l-L+(/3)/L+(fc) 
+ iA; 2 V + i{k-P) 



+ 



l-L+(-/3)/L+(fc) 



+ i{k + f3) 



) 




+ 



+ i{k + l3)) 



G{k) . (62) 



In particular, it follows that 



^(+0) = w{-0) = w{0) = Co L+(0) - ^ + L+{-/3)) 



(63) 



We stress that the solution corresponds to the limit a — 0. 



4 The critical displacement w{0) 



We present in this section the dependance of the critical displacement w{0) on the velocity of 
failure propagation v. This gives precise information on the range of velocity where steady- 
state regime is attainable or not, and where steady-state propagation is stable or unstable, 
information that cannot be obtained from the uniform continuous model of Sect.O Also, 
it is possible to establish a connection with subsonic, intersonic and supersonic regimes of 
the uniform continuous model considered in Sect. [2] and find a universal character of the 
propagation which does not depend on some geometric parameters as the pillar normalized 
distance a. In this work we focus our attention on numerical results where ml = 
Ml = M2 and X2/X1 < 1; the analysis of the model, where a contrast between the inertial 
properties of the damaged and undamaged structure is present, is left for a future work. 

In the limit a = +0, ArgL(A;) is a piecewise constant function over a semi-infinite support 
{v > 0), which makes the computation of the integral on an unbounded domain in Eq. fl60l) 
particularly difficult. To overcome this problem we consider a simplified model with massless 
beams where the inertial effects are concentrated in correspondence of the elastic supports 
(shown in Fig. [lb) and a second one where we introduce a small positive dissipative term a 
(see Eq. (HH])). 

4.1 The lattice beam model for masses concentrated at nodal 



We evaluate the displacement according to formula ( !59|) . which also include evaluation of 
the Arg[L(/c)] and of the integral in Eq. ( I6OI) . This results will enable us to judge on the 
values of the speed of the propagating fault in the discrete flexural system. We consider the 
case when the normalised masses Mi and M2 used in the previous section are now replaced 
by the total redistributed mass, which includes the added masses assigned to junction points 
X = na and the mass assigned earlier to the beam between the neighbouring supports, as (in 
the non-dimensional form) Mi+m^ /rrii = 1 and M2+m°/mi = 1 — M1 + M2, respectively. In 
this section, we assume that Mi and M2 are equal, and we use the notation M = Mi = M2. 



points 
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As a result, we deal with a system of unit point masses placed at x 
non-inertial elastic beams. In this case the equation f p8|) reads as 

d'^w{x, rj) 



an and connected by 



0. 



and hence, when an < x < a{n — l), the displacement is cubic in x. For x 
are the same as in (l36l) and (!37l) . 



(64) 

an the equations 



UJ 



VlV^^Viii Vi 




Figure 10: Dispersion diagrams uJi{ka), x = xi, and U2{ka), x = X2, for waves in the 
discrete flexural system. Fig. [lb. The diagrams are given for stifness ratio K2I >c\ = 0.5 
and normalized span length a = 2.0. The straight line u = kv {v = 0.23) has a finite 
number of intersections (denoted with small black circles) at wavenumbers k^\ . . . , k^^ and 
k^\ . . . , k'^^ with the two dispersion curves uJi and UI2, respectively. The velocities vj, vjj, 
Viii, Vjv are reported in the critical displacement diagram in Fig. 



The Wiener-Hopf equation for the one-sided Fourier transform of the displacement has 
the form (HSil with Mi = M2. Taking into account Eq. (jHj) we deduce 



Li{k) = 1 + {0 + ikvf + 



2 ^ 12 (1 — cosfca)^ 
^ 2 + cos ka 



L2{k) = X2/X1 + {0 + ikvf 



12 (1 — cos ka) 



a^ 2 + cos ka 

Consequently, the dispersion equation describing the waves for 77 > and 77 < are 



(65) 



12(1 - cosA;a)2 
a^(2 + cos ka) 
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X2 12(1 — COS kay 
^^1^1^ a4(2 + cos ka) 



(66) 



respectively. The corresponding dispersion curves are shown in Fig. [TOl 
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Figure 11: The piecewise constant (a) Arg[Li], (b) Arg[L2] and (c) Arg[L] as a function of 
the normahzed wavenumber ka. Results are given for a = 2.0 and v = 0.23. Jumps occur at 
the intersection point wavenumbers k[^^ , • • • , ^3^"* and k^^^ , • • • , ^5^^ highlighted in Fig. [TUl 

Due to Eq. (!52l) LQ{k) = L{k) = Li{k)/L2{k). In this case, Aig[L{k)] is a piecewise con- 
stant function with a finite, speed- dependent support {v > 0) vanishing in the neighborhood 
of the origin. The jumps occur at the points of intersection at A; = kj^\ k = k^p between the 
ray u = kv and the dispersion curves uji^2{k), as shown in Fig. [Till The graphs of Arg[Li(A;)], 
Arg[L2 (/;;)] and Arg[L(/i;)] for v = 0.23 are given in Fig. [11] According to Eqs. ([59]) and (!60|) 
the displacement at the origin is given by the formula 



w{0) = g 



2n.2+l 




'^2ni+l 



(67) 



where 2ni + 1 and 2^2 + 1 are the number of intersections of the straight line u = kv, with 
the two dispersion curves u = uJi{k) and u = u}2{k), respectively. Note that in Fig. [TU] and 
Fig. [m at = 0.23, Ui = 1 and na = 2. 

The critical displacement w{0) as a function of the speed v is given in Fig. [12] for the 
stiffness ratio xa/xi = 0.5 and span length a = 2 and compared with the continuous beam 
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w{0) 




Figure 12: Critical displacement w{0) as a function of the speed v. The uniform continuous 
model (Fig. [T^, eqs. [181 ED and [25]) is compared with the discrete massless beam model 
(Fig. eq. 1^7]) . Results are given for K^j ycx = 0.5 and in the discrete model a = 2. The 
velocities vi, vu, vjjj, vjv are reported in the dispersion diagram in Fig. [101 The critical 
displacement Wc/g = 0.3 is denoted by the grey line. 

model on elastic foundation. It is evident the jagged structure for v < vi, where the straight 
line u = kv has multiple intersection with the dispersion curves ui and/or U2 and, in 
principle, different waves can radiate from 77 = 0. At w = 0.23 waves can radiate with 

(1) (2) (2) (2) 

wavelength = in the region rj > and with wavelengths = fc^ , fcg , fc^ in the region 
?7 < (see Fig. [TOl) . In this low-speed region more and more sinusoidal waves arise as the 
speed decreases, with intensities extremely sensitive to the speed, causing irregularities in 
the dependence of the critical displacement on v. 

However, additional conditions have to be fulfilled by displacement fl^Tl) to be a steady- 
state solution of the problem f[5B]) - f[571) . As noted by Marder and Gross (1995) for the problem 
of dynamic crack propagation in lattices, the solution is linearly unstable whenever w{0) is 
an increasing function of v, which means that the failure wave accelerates. Therefore, all the 
increasing portions of the curve can be ruled out. 

But this is not enough, a further condition of admissibility needs to be satisfied by w{0): 
it is not only necessary to verify that the displacement w{0) has reached the critical value 
Wc (positive in normalized coordinates) at t = 0, but this must be the first time at which 
the displacement reaches this critical value. This is equivalent to check that w{r]) < Wc for 
?7 > 0. Instead of looking at the distribution of the displacement for 77 > 0, Slepayn and 
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Ayzenberg-Stepanenko (2004) considered a necessary condition equivalent to 

dwiri) 1 dw , , , , 

— ^ = — — <0 r/ = 0, 68 
drj V at 

for the problem of transition waves in bistable-bond lattices. Considering the limiting relation 

w{0) 



w{0) = lim (— ifc)(i/cf) 

fc— >ioo 



(69) 



and that for k — )■ +ioo 

1 1 1 P lni^(0 



L+{k) 2niJ_^ i-k 

we check the range of velocity where 



l-TT^/ ^-^de, (70) 



u;(0) = / lnL(A;) dA; > 0. (71) 



2tx 

By using the symmetries of Li{k) and L2{k), the integral in (17T]) is given by 

/OO /'OO 
\nL{k)dk = 2 \n{\L{k)\) dk. (72) 
-OO Jo 

which is integrable but the integrand has logarithmic singularities. Then, the integral can 
be split as 

/ + / (73) 

where (/ci, /c2, . . . , /c2(ni+n2+i)) are the 2ni+2n2+2 intersection points between the line u = kv 
and the dispersion curves cui and U2, such that ki < kj for i < j. For the parameters con- 
sidered in Fig. [ini at v = 0.23, the wavenumbers are {k^\ kf \ kf \ k^^\ k!^\ kf\ kf \ k^^^). 
The first 2ni + 2n2 + 2 integrals can be computed numerically using standard quadrature 
rules for integrand having logarithmic singularities at the endpoints and /cm has been chosen 
to be sufficiently large to guarantee that | \n{\L{k)\) dk\/\ J^"^' ln(|L(fc)|) dk\ < 10"^. 

These values are reported in Table [H We see that for v < 0.55, w{0) < (except at 
V = 0.2), showing that in this range no steady-state propagation, w = w{x,ri), can exist. 

Following the behavior of the transversal displacement w{0) at higher velocities, it can 
be seen that there is a sort of plateau in the range vj < v < vn with a drastic drop in the 
range vn < v < vm. This dependence shows strong similarities with the behavior of the 
continuous model. 

We can conclude that steady- state propagation, when possible, occurs in a range of velocity 
which is exceptionally well approximated by the narrow range of velocities ^MEj, within the 
intersonic regime, provided by the continuous model of beam on elastic foundation. In such 
a velocity range, the steady state solution for the discrete case is linearly stable and satisfies 
the conditions (17T!) and the critical one w(0) = Wc- 

In Fig. [13] the displacement w(0) is given as a function of velocity v for different values 
of the normalised span length a. We note the good agreement between the continuous and 
the discrete models for sufficiently small values of a, which is expected on physical ground, 
an issue that will be better discussed in Sect. 14.21 Also, the continuous model is approached 
from below at decreasing values of a. 
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0.31 * 10-2 




1.00 


0.52* 10-1 


0.25 


-0.12* 10 
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0.62* 10-1 


0.30 


-0.80* 10 
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0.24* 10-1 
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0.77* 10-1 


0.45 


-0.18* 10 


-1 




0.85 


0.30 * 10-1 




1.2 


0.11 * 10-1 



Table 1: Value of the displacement time derivative w;(0) at different speeds v. 




Figure 13: Critical transversal displacement w{Q) as a function of the speed v. The discrete 
massless beam model (continuous grey lines) and the uniform continuous model (dashed black 
line) are represented with X2/X1 = 0.5. In the discrete model a = 0.25, 1.0, 2.0, 4.0, 8.0. 

4.1.1 Dispersion diagram transformation and the continuous limit 

Let a — > +0. Expanding in ka we deduce from (166|1 that in any finite range of the wavenumber 
the dispersion relations take the form 

~ Vl + fc4 , UJ2 ~ Vxs/xi + , (74) 

while the remaining periodic part of the dispersion curves moves away to infinity (see Fig. fT^ . 
In this transformation, the two branches, Ui and U2 fl66|) . coalesce, becoming closer and 
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Figure 14: Dispersion diagrams u!i^2{k) for the discrete massless beam model. Results are 
given for a = 1.4, 1.2, 1.0, +0 and K^j hx = 0.5. 

closer, and the separation uJi{k) — U2{k) as a 0. Hence in the limit the cross points 
other than those corresponding to the asymptotic formulae ( 17^ . do not contribute to the 
product in ( 1671) for the displacement w{0). At the same time, these asymptotic expressions 
coincide with the dispersion relations (fT5|l derived for the continuous model, which has the 
same average density and the same average support stiffness as the discrete model. Thus, 
the dynamic response of the discrete fiexural system is shown to converge to the continuum 
limit as a — )■ 0. 

4.2 Inertial beam model 

For the case of dynamic beam on the set of discrete elastic supports (Fig. [lb) the dispersion 
diagrams for the undamaged and damaged structures have infinite numbers of curves (see, 
for example, Brun et al., 2012) and Eq. ( 1671) is still valid, but with ni, n2 — )■ oo. To compute 
the displacement w{0) is then convenient to introduce some dissipation as in (149|) and make 
use of Eqs. (1591) and (|60I). 

In Fig. [15] the displacement w{0) is given as a function of velocity v for different values 
of the normalised span length a. Computations have been done assuming a = 0.01 and the 
integral in ( l60l) was truncated imposing a relative error smaller than 10~^. The comparative 
analysis between the results of Fig. [T3]and Fig. [T5] shows larger differences at increasing values 
of a, where the difference between the coefficients in ( H3|) and in (jH]) is more pronounced 
and the wave propagation is better described by the beam inertial model. 

Despite of the fact that the interval of velocity where oscillations are present increases 
with a, all the curves display a drastic drop in correspondence of the intersonic regime 
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Figure 15: Critical displacement w{0) as a function of the speed v. The discrete inertial 
beam model (continuous grey lines) and the uniform continuous model (dashed black line) 
are represented with X2/X1 = 0.5. In the discrete-continuous model m^/rrii = 0.67 and 
a = 0.25, 1.0, 2.0, 4.0, 8.0. 

identified by the model of beam on elastic foundation, thus giving a universal property for 
the steady-state failure propagation in such mono- dimensional structures. 

Finally, we report some data deduced from real life bridges, showing the range of possible 
engineering applications of the proposed model. We first consider a bridge of average dimen- 
sion, the S'Adde bridge, whose geometrical and material properties are given in the caption 
of Fig. 4 of Brun et al., 2012, with a span length of 90 m. In such a case the normalised 
span length is a = 2.02 for vertical flexural waves and a = 0.83 for horizontal flexural waves. 
In addition, we also pay attention to the Millau viaduct, one of the tallest vehicular bridge 
and the longest multiple cable-stayed bridge in the world (Magalhaes et al., 2012). The steel 
structure has Young Modulus E = 210000 MPa, moments of inertia Jy = 1137.7 m^ and 
Jz = 8267.5 m^, vertical and transverse stiffness = 20000 MPa m and >iy = 41.84 MPa 
m, respectively, and span length of 342 m. Then, the structure has normalised span length 
a = 7.61 for vertical flexural waves and a = 0.99 for the horizontal ones, showing that, also 
for this extreme structure, we are in the range of a where the presented model can be used 
to model the steady-state failure propagation. 
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5 Concluding remarks 



We have obtained novel results on an advance of a transition flexural wave through a beam- 
like periodically supported discrete structure. Both propagating and evanescent waves are 
included in a general solution and hence influence the interaction between different nodal 
points within the discrete system. 

It was shown that in the considered problem, in contrast to the 'classical' fracture in 
continuum mechanics, the energy release under the support damage remains positive not 
only in the case where there is no wave radiation but also in a part of the intersonic speed 
range. In this case, a part of the total energy released from the heavy beam is radiated by 
an elastic wave excited by the jump in the stiffness of the support. 

The model has several significant applications in structural mechanics, including safety 
evaluation and design algorithms for earthquake resistant long bridges. The necessary condi- 
tion, reduced to evaluation of the critical displacement, under which a failure of the support- 
ing piles of a bridge can propagate, has been identified and analysed for a range of physical 
parameters. The failure has been modeled as a drop of the stiffness or of the total mass 
the support. For the model of inertial continuous beam on distributed elastic foundation 
the solution has been obtained by solving directly the corresponding equations of motion 
for steady-state propagation and checking, in addition, the energy balance. Computational 
examples, as outlined in the above section, strongly suggest the applicability of the proposed 
model to real life bridge systems such as the S'Adde bridge and the Millau viaduct. 

We note that the model discussed in the paper is radically different from models describing 
crack propagation in discrete lattices, both in scalar (anti-plane shear) and vector cases. The 
present models deal with the fourth-order flexural wave equation, whose dispersion properties 
include unique features, which are absent in the two-dimensional lattice crack models. 

A regime, which we refer to as intersonic, has been identified (\/2(x)^''^ < v < V2), and 
it has been shown that the solution corresponding to the steady propagation of the interface 
wave is intersonic. An efficient and elegant mathematical approach based on analysis of a 
functional equation of the Wiener-Hopf type, has led to the expression of the displacement on 
the interface wave to be written as a product of terms evaluated directly from analysis of the 
kernel function of the Wiener-Hopf equation. The discrete model has been compared with 
the simpler continuous approach, and of course it has been shown that the lattice solution 
possesses new feature, not reported in the past for interfacial flexural waves. 

Lastly, the conditions allowing the failure wave to propagate and the characteristic wave 
speeds are found. 
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